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Abstract 

We extend the usual chiral perturbation theory framework (xPT) to allow the inclusion of 
a light dynamical isosinglet scalar. Using lattice QCD results, and a few phenomenological 
inputs, we explore the parameter space of the effective theory. We discuss the S-wave 
pion-pion scattering lengths, extract the average value of the two light quark masses and 
evaluate the impact of the dynamical singlet field in the low-energy constants /i, /s and ^4 of 
xPT. We also show how to extract the mass and width of the sigma resonance from chiral 
extrapolations of lattice QCD data. 
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1. Introduction 



Chiral Perturbation Theory, xPT [1,2], has become a standard tool for the phenomenological 
description of QCD processes involving pseudo-Goldstone bosons at low-energy (see [3] for 
a review). It is grounded in a few simple assumptions: (i) the underlying theory of strong 
interactions, namely QCD, has an exact SU{nf) x SU{nf) chiral symmetry in the limit of 
vanishing Uf light quark masses that is spontaneously broken down to SU{nf), (ii) there is 
a mass gap (~ A^) for all states except for the Goldstone bosons, and (iii) the exact chiral 
symmetry is explicitly broken by the actual non-vanishing quark masses, rriq <^ A^. Under 
those assumptions one can construct a low-energy effective theory, xPT, for the Uf pseudo- 
Goldstone bosons organized in powers of p/A^ and m^/A^, where p is the typical momentum 
of the low-energy processes {p <^ A^). In practice, A^ is taken of the order of the rho mass 
{nip ~ 770 MeV), and the pseudo-Goldstone bosons are identified with the pions for n/ = 2 
(?Ti7r ~ 140 MeV) and with the lightest octet of pseudoscalar mesons for Uf = 3, which also 
includes the kaons {rriK ~ 490 MeV) and the eta (m^ ~ 550 MeV). 

Scattering amplitudes can be systematically calculated within this framework to a given 
order in ~ over A^ . However, when pion scattering amplitudes are calculated in the 
isoscalar channel, a bad convergence is observed, even at reasonably low-momenta. This 
has led some authors to resum certain classes of diagrams, using a number of unitarization 
techniques (see, for instance, [4, 5, 6, 7, 8]). Most of these approaches improve considerably 
the description of data with respect to standard xPT, and indicate that a scalar isospin zero 
resonance at relatively low-mass, the sigma, exist. In fact the mass and width of the sigma 
resonance are nowadays claimed to be known very accurately rria = 44ll8^ MeV , r/2 = 
272+^2,5 MeV [9, 10] (see also [11]). 

Under the SU (3) perspective one may find surprising that the effective theory contains 
kaons but not other states with similar masses, but different quantum numbers, that can be 
equally excited in a collision at intermediate stages. The relatively low-mass of the sigma 
resonance, with respect to the chiral cutoff, A-^, and its proximity to the value of the kaon 
mass suggests that it may be convenient to introduce it as an explicit degree of freedom in an 
extension of xPT, thus lifting the assumption (ii) above. It is in fact an old observation by 
Weinberg [12], that the explicit inclusion of resonances in a Lagrangian generically improves 
perturbation theory. 

We implement this observation here in a chiral effective theory framework that involves 
a dynamical singlet field together with the lowest pseudo-Goldstone bosons. We write down 
the most general SU (2) chiral Lagrangian including an isospin zero scalar field at order 
and calculate a number of observables at this order. We show that for a large scalar mass, 
the effect of the scalar reduces to just redefinitions of the low-energy constants (LEG), hence 
explicitly demonstrating that our approach is compatible with standard xPT. However if we 
count the mass of the scalar as order p^, namely of the same size as the pion mass, the non- 
analytic pieces of our amplitudes differ from those of xPT. Furthermore, the quark mass 
dependence of the observables is also different. We compare this effective theory, which we 
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call xPT^, versus standard xPT against lattice data on the pion mass and the pion 
decay constant F^^ [13], and on the pion-pion S-wave scattering lengths [14]. At the current 
precision the lattice data is unable to tell apart xPT from xPT^. 

We organize the paper as follows. In the next section we discuss the power counting, 
construct the Lagrangian up to next-to-leading order (NLO), and compare it with the one 
of the linear sigma model. In section 3 and section 4 we calculate the two-point function 
of the axial- vector current and of the scalar field respectively, up to NLO. In section 5 we 
perform a number of fits to lattice data for m^r and both in xPT and in xPTs in order to 
constrain the parameter space of the latter. In section 6 we discuss the LO S-wave scattering 
lengths, and compare the results of xPTs together with those of xPT with very recent lattice 
data. In particular it is shown how the mass and decay width of the sigma resonance can 
be extracted from them. We close with a discussion of our results and the conclusions in 
section 7. 



2. Lagrangian and power counting 

Our aim is to construct an effective field theory containing pions and a singlet scalar field 
as a degrees of freedom, that holds for processes involving only low-energy pions as the 
asymptotic states 

p , m^(~ 140 MeV) , ms{-^ 440 MeV) < A^. (2.1) 

The structure of the effective Lagrangian will be independent of the underlying mechanism 
of spontaneous chiral symmetry breaking. It consists of an infinite tower of chiral invariant 
monomials combining pions and a singlet scalar field with the generic appearance 

= ^ ^{k,l,r) , (2.2) 
(k,l,r) 

where jC(^k,i,r) contains k powers of derivatives, / powers of the scalar or pseudoscalar sources 
and finally r powers of the singlet field. 

being p a typical momentum in the process that we have assumed to be of the order of 
the scalar particle mass. One possible manner to relate these scales is to assume that 
~ ml ~ ''TiqA^, like in standard xPT- Hence, in the chiral counting C[k,i,r) is of order 
pk+r+2i^ In fact, in this paper we only use the inequalities in (2.1). More refined hierarchies, 
like <^ , p <^ A^ may be interesting to explore in the future. Notice that terms with 
k + r + l < 4 correspond to relevant operators and, hence, their dimensionful constant may be 
tuned to a scale smaller than the natural one A^, as it happens in standard xPT {F^^ <^ A^). 
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The Lagrangian involving pions and scalar fields transforming as a singlet under 5'f/(2)/j x 
SU(2)l, respecting Chiral symmetry, P and C invariance, has been presented in the linear 
approximation in [15, 16] and up to quadratic terms in [17]. For the time being, we will collect 
only the relevant terms necessary for our purposes. The leading order (LO) Lagrangian 
consist of three parts: the standard Goldstone boson Chiral Lagrangian, that we do not 
discuss further, terms involving the scalar field only, and interaction terms between the 
scalar field and the pseudo- Goldstone bosons. 

2.1. Leading Lagrangian 

Consider first the part of C^^ containing only the singlet scalar field. In the absence of any 
symmetry hint we are forced to write the most general polynomial functional, 

= ^d.Sd^S - ^m'sSS - XiS -^S'-^S' + --- (2.4) 

where the dots indicate terms suppressed by powers of 1/A^. Suppose that we deal with the 
chiral limit. At LO Ai must be set to zero in order to avoid mixing of S with the vacuum, 
and at higher orders it must be adjusted for the same purpose. The mass and the coupling 
constants above are functions of the small scale ths and the large scale A^, (m^ <^ A^). Their 
natural values would be A3 ~ ^{^x) ^^"^ -^4 ~ ^(1)- ^^at case, the scalar sector above 
becomes strongly coupled. However, strongly coupled scalar theories in four dimensions are 
believed to be trivial [18, 19]. Their exact correlation functions factorize according to Wick's 
theorem and consequently they behave as if the theory were non-interacting. A practical 
way of taking this fact into account is just setting A3 = A4 = 0, which we will do in the 
following. When the interactions of the scalar with the pseudo-Goldstone bosons are taken 
into account, small (m|/A^ suppressed) but non- vanishing values of A3 and A4 are required 
to ensure perturbative renormalization of the whole C^^. 

The second contribution we are interested in is the lowest order Lagrangian describing 
the interaction of the scalar field with the pseudo-Goldstone bosons. As a basic building 
block we use the unitary matrix U (x) to parameterize the Goldstone boson fields, that may 
be taken as, 

although final results for observable quantities do not depend on this specific choice. At LO, 
F may be identified with the pion decay constant F^. We also use the building block, 

X = s + ip = 2Bml . (2.6) 

On the r.h.s. we have set the pseudoscalar source, p, equal to zero and s to the diagonal 
matrix. As we will work in the isospin limit we will use rh referring to the average quark mass 
between ruu and m^. The covariant derivative acting on U (x) is defined as usual, containing 
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the external vector and axial sources, D^U = dfj^U — i[v^, U] +i {a^, U} . The transformation 
laws for all these building blocks under the local symmetry group G = SU{2)j^ x SU{2)l 
are dictated by 



Qr gR^ + ignd^gR^ , = - gil^^gL^ + igLd^g^^ 



(2.7) 



where gR, gi C SU (2) . 

With all those ingredients one can construct the Lagrangian, 

+ (t + ^""'"^^ + + ■ ■ ■ ) ^^^^ + ^^^^ ^^-^^ 

where the ellipsis stand for higher order terms involving more powers of the singlet field (or 
derivatives on them), which are suppressed by powers of 1/A^. 

At this point a small digression is in order; notice the peculiarity of (2.8) with respect to 
the usual chiral expansion. At this order, both expansions can be cast in the form, 

>C-^E^'^(^x,5)0« (2.9) 

k 

being an operator of order k including only the pseudo-Goldstone bosons and bk{A^, S) 
its corresponding "Wilson coefficient" , that can depend on the singlet field if one considers 
the theory with the scalar field inclusion. While in the standard theory the power counting 
is given entirely by the operator, i.e. bk{A^) ~ (9(A^~'^), in the extended version one also 
has to take into account that the Wilson coefficients themselves have a power expansion in 
S/Ay.. At higher orders operators containing the derivatives of the scalar field must also be 
included. 

Before closing this section we would like to remark that even if we have kept for F and B 
the same names as in xPT, they are now parameters of a different theory and, hence, their 
values are expected to differ from those in xPT. 



2.2. Comparison with the linear sigma model 

Hitherto we have included in a dynamical fashion a scalar particle interacting with pseudo- 
Goldstone bosons. One may wonder if there is any relation between the effective theory just 
introduced and the old linear sigma model [20] (see [21] for a review), which we discuss next. 
The starting point for the construction of the linear-o" model is an 0(4) invariant action. 
The global 0(4) symmetry is spontaneously broken down to 0(3) because the scalar field 
develops a non-zero vacuum expectation value v. 
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The Lagrangian reads 

= ^ {d.ad^a + d,^d^0) - ^ + _ ^ (2.10) 

where (p is an isotriplet pseudoscalar field, usually identified with the pion, and a is an 
isosinglet scalar field that, after the shift a — )■ a + f is usually identified with the sigma 
resonance. We do not display the part of the model containing nucleons because it has no 
relevance for our discussion. Since 0(4) = SU{2) x SU{2) for group elements close to the 
identity, the model transforms correctly under the SU (2) ^ x SU {2)l chiral symmetry of two- 
fiavor QCD with massless quarks. To see this explicitly we make the change = all — ir ■ 0, 
being r the Pauli matrices. Then (2.10) can be written as 

= hd.^d^J:^) - ^ ((StS) - 2v'Y , (2.11) 
4 lo 

that explicitly exhibits the desired symmetry (2.7), if we transform gj^T, g^^ [3]. The 
traditional identification of the ip fields with the pions and the a field (after the shift) 
with the sigma resonance, which is fine concerning the transformations under the unbroken 
subgroup 0(3) = SU{2), becomes problematic if one wishes to implement the non-linear 
SU{2)ii X SU{2)l symmetry that the model retains after the shift a — )■ cr + w is performed. 
In order to make the non-linear SU{2)ii x SU{2)l symmetry manifest in the Lagrangian 
above and keep the transformations of the Goldstone bosons in the standard way [22, 23], 
as we have done in the previous section, it is convenient to perform a polar decomposition 
of S, E = (t" + S)U , with U being a unitary matrix collecting the phases, to be identified 
with the U appearing in (2.5), and S a real scalar field, to be identified with our singlet field 
above. We remark that S must not be mistaken by the a field in the original variables of 
the linear sigma model. The symmetry transformations of the fields S and U are the same 
as in (2.7). This change of variables leads to 

= + |5 + Is'^ {D,UD>^U^) + Id.Sd-S - Xy [s' + ^ + ^)- (2-12) 

The terms with covariant derivatives above have the very same functional form as the terms 
with derivatives of (2.8), with the identifications v = F, cm = 1/2 and C2d = 1/4. However, 
the terms with no derivatives, the potential, are set to zero (or, at higher orders, to small 
values uncorrelated to the rest of the parameters) in xPT^, except for the mass term, for 
which m| = 2Xy. This is because the underlying mechanism of chiral symmetry breaking 
is assumed to take place at the scale A^, and hence it must not be described in the effective 
theory. 

Since pions are not massless in nature, a small explicit breaking of the 0(4) symmetry 
had to be introduced. This was traditionally done by adding a term SC^ = Ha. In terms of 
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the new variables this term reads 

5C^ = Ha = + S^) = ^{v + S){U + f/t) . (2.13) 

Hence, it has exactly the same functional form as the terms with no derivatives in (2.8), 
once X is set to 2i?mll, with the identifications H = 2FBm, = 1/4 and C2m = 0. 

In summary, the Lagrangian of xPT^ at LO differs from the one of the linear sigma 
model only in two respects: (i) the self-interactions of the scalar field S are set to zero (or, 
at higher orders, to small values uncorrelated to the rest of the parameters), and (ii) it has 
four additional free parameters controlling the interaction of the scalar field S with the pions: 

Cld, C2d, Cimi C2m- 



2.3. Chiral symmetry constraints 

To envisage the effects of explicit chiral symmetry breaking on the dynamics of the singlet 
field we set U to the vacuum configuration {U = 11). The terms proportional to the quark 
masses in (2.8) induce new terms in the Lagrangian of S*, that can be reshuffled into the 
coefficients of (2.4). For the first two terms one finds explicitly 

Ai — )■ Ai — SFcimBm , m| — m| = m| — lQc2mBrh . (2-14) 

As a consequence the singlet field is brought out of its minimum in the chiral limit by terms 
proportional to m. Hence, the direct consequence of the inclusion of non-vanishing quark 
masses results in a new contribution to the singlet-vacuum mixing. The new scalar field 
describing the first excitation with respect to the vacuum may be obtained by carrying out 
the following shift 

S^S + FSo with 5o = 8ci™^-^. (2.15) 

rrig mi;}' 

After this shift, and upon separating the vacuum contribution, the original Lagrangian (2.8) 
keeps essentially the same form, 

= (^^rod + Fr.dS + r2dS' + ■ ■ ■ ^ {D^UD^U^) + 

^r^m + Fn^S + r^n^S^ + ■ ■ {{x^U + xU^) - (x^ + x)) , (2-16) 

provided we redefine the LEG as 

rod = 1 + 4cidS'o + 402^5"^ + . . . , ro„ = 1 + 4ci„S'o + 4c2m,S'o + • • • , 

rid = Old + 2c2dSo + . . . , ri„ = ci„ + 2c2m5'o + . . . , (2.17) 

r2d = C2d + • • • , r2m = C2m + ■ ■ ■ ■ 
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In the previous expression all the terms explicitly depicted are C(l) quantities and ellipsis 
stand for subleading contributions CnxSo"'~^^ ~ (F/A^)"~^, for n > 2 (x = d,m). 

There is a subtle point that must be addressed before going on: for generic values of the 
LECs the shift (2.15) breaks chiral symmetry. This is most apparent if we lift the scalar and 
pseudoscalar sources from its vacuum values to arbitrary ones. Namely, if the original scalar 
field in (2.4) is a singlet under chiral symmetry, the scalar field after the shift (2.15) is not. 
This is so for any value of the parameters, except for those that fulfill 

Ai = ^i5^. (2.18) 

If we choose Ai as above, the shift becomes independent of the quark masses (5*0 = — Cim/2c2m) 
and hence the scalar field after the shift is still a scalar under chiral symmetry, as it should. 
Moreover, for this choice, rim = 0, and tq^ and rod can be set to 1 by a redefinition of B 
and F respectively. The net result is equivalent to choosing Ai = cim = in (2.4) and (2.8). 
This has in fact a simple interpretation. If we impose to our original scalar field to be a 
singlet under chiral symmetry for any value of the external sources and not mix with the 
vacuum, then the only solution at tree level is Ai = Cim = 0. We shall adopt this option 
from now on. At higher orders these two parameters must be tuned so that no mixing with 
the vacuum occurs at any given order. Note finally that the value cim = is incompatible 
with the linear sigma model one cim = 1/4. 



2.4. Next-to-leading Lagrangian 

In computing loop graphs, we will encounter ultraviolet divergences. These will be regular- 
ized within the same dimensional regularization scheme as used in [2], and the elimination 
of the divergences proceeds through suitable counter-terms. Like in [2] we will deal with 
contributions up to including terms of order p^. Since the singlet fields will only enter in 
internal propagators, the counter-term Lagrangian we need only involves pions, and hence 
it has the same functional form as the one in xPT. The coefficients, however, receive extra 
contributions due to the appearance of the Ci^ {i = 1,2, x = d,m) bare parameters. In 
addition, unlike standard xPT, now F and B need to be renormalized. In order to take this 
into account we chose to include explicitly the corresponding counter-terms below. Using 
the SUl{2) X SUr{2) formalism [24], rather than the 0(4) one [2], we have 

£(4) = ^iiiD^UDf^U^y + U2{D,UD^'U^){D''UD^U^) 

+ ^HX'U + xU^? + \h{D^U^D^x + D'^X^D^U) (2.19) 
+Ziml{x^U + xU^) + Z^ThliD^UD^'U^) + .... 
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In order to avoid confusion with the values that the LEG take in xPT and xPTg, we shall 
denote the former /j, -i = 1, . . . , 4 and the latter £j, i = 1, . . . , 4. The relations between /j, /[ 
and li, i = 1, . . . ,4 that appear in the paper are the standard ones in xPT [2]. In this work 
ii and £2 will not be necessary for renormalization. For the observables that we will consider, 
the pole at = 4 is removed by the following two kinds of renormalization constants which 
occur in the Lagrangian C^^^ 

ir-ei + 'lA (i = 3,4), Z,. :=ZJ + r,A = 1,2), (2.20) 

The first ones, 7^, are a simple redefinition of the divergent part in the standard monomials 
of xPT 

73 = -^ + 32C2m(c2m-C2d)-8Ci^(l -4C2m) , 74 = 2 + 40^^ (1 - 8C2m) + 32C2dC2m • (2.21) 

While the second, absent in xPT, are entirely due to the interaction of pions with the singlet 
field. In our case, this contribution has two sources; one proportional to B and other to F 

Ti = -2 {cj^ - C2d + C2m) , r2 = C?, - C2d ■ (2.22) 

Recall that, like F and B, and £4 are now parameters of a theory different from xPT 
and, hence, their values are expected to differ from the ones in the latter. Note also that 
our approach differs of that in [15] in the respect that pions and the singlet field are both 
dynamical in the same energy range and thus both will be allowed to run inside loops. In 
that respect, any estimate of the LEG by matching the observables derived from xPTs with 
those obtained from a Lagrangian of resonance exchange as in [15] should keep the singlet 
field S in the latter as a dynamical low-energy degree of freedom. Loop effects of scalar 
resonances coupled to pseudo-Goldstone bosons have been studied in [17, 25, 26, 27]. 



3. The axial-vector two-point function 

We are now in the position to perform a complete NLO analysis of the pion mass and decay 
constant, including the radiative correction due to the singlet field. In order to calculate 
them we focus in the sequel in the axial-vector two-point Green function. The first quantity 
will be defined through the position of the pole in the two-point function while the second 
one can be obtained directly from its residue. 

At the Born level the expressions for the pion mass and decay constant do not differ from 
those of the standard xPT theory while at NLO corrections are slightly more cumbersome. 
We denote by mpg and Fps the pion mass and the pion decay constant respectively calculated 
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(a) (b) (c) (d) 

Figure 1: Self energy diagrams. Solid lines denote the scalar and the dotted ones pions. 
Diagram (a) corresponds to NLO counterterms 




(a) (b) (c) (d) 

Figure 2: Diagrams contributing to the decay constant. Solid lines denote the scalar and 
the dotted ones pions. Diagram (a) corresponds to NLO counterterms. 



at NLO, whereas we keep m,r = 2i?m and F for the same quantities at LO. The diagrams 
contributing at 0{p'^) to mpg and Fpg are represented in Fig. 1 and in Fig. 2 respectively. 
In both figures the diagram (a) is the usual couterterm contribution, and the diagram (b) 
the usual tadpole contribution, already encounter in the standard theory. Diagrams (c) and 
(d) are new and appear because of the singlet field. 

We cast the expressions for the mass and decay constant as 

mis =25m + Um. + Pm + 0{p') , 
Fps=F{1 + Uf + Pf + 0{p'')) . 

They contain contributions related to the unitarity cut {Um, Up) in the s-channel. Fig. 1(d) 
and Fig. 2(d), and a polynomial term in s (= m^) which includes logarithms {Pm, Pp)- 

Notice that, as we introduced a new degree of freedom, the infra-red physics is modified 
and in particular the analytic structure of the amplitudes. While xPT finds its first cut at 
NNLO for the pseudoscalar two-point function, located at 3m7r, xPT^ has already at NLO 
a cut located at ms + This enhancement can be understood if one chooses to interpret 
the scalar field as a simulation of a strong two-pion rescattering. Consider for example the 
three pion contribution to the pion self-energy. If we take the two-pion subdiagram, say 
t(^), and replace it according to t^^^ (t(2))2/ (^(2) _ ^(4)j _ ^(2)^ 

as suggested by the Inverse 

Amplitude Method (lAM), the first term would correspond to the pion-sigma contribution 
in XPT5 and the second one to a pion tadpole, both contributions being of O (p^) in the 
pion self-energy in xPTg. 
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The expressions for the mass are given by 



4^2 



Ami f fJ'S - l^n 



and for the decay constant 



2r2 



— J {ml, ml; ml) {m^ - 2mlY , 



Cld^S - 4:C2m^l^ns) + 



4 2 ° 2 

ml - mtmi; 



(3.2) 



1 1 2, I A 

J{m^,ms;m^ 



Ami -ml I ^ 



3m|) , 



?j (ml — 2m'i)^ Ami 



^id 



g^2p2 _ j^2 



■> 2 

m 



;^7r - /i5 

r2.^2 • 



clAml-2mlf 

—Ta 5 2\ 1" ^C2mi 2^TTi 

[Amf^ — mg) 



(3.3) 



The functions J and /i^ {a = n , S) are displayed in the Appendix. In addition we have 
used the scale independent quantities ii and Zj, that are defined as follows: 



7i 



327r2 



3,4) 



327r2 



(J = l,2), (3.4) 



where A is the renormalization scale. 

Both quantities in (3.1) have the following virtues that constitute non-trivial tests on 
their correctness: 

1. Despite their appearance, they are finite in the chiral limit, m — t- 0. More explicitly, 
in this limit the pion mass vanishes, as it should, while the decay constant reads 



PS 



F 1 + 



Z2T2 + -Cirf 



(3.5) 



2. Setting — )■ (x = m, (i) in (3.1) they reduce to their standard xPT values 



^ps xPT = ( 1 - Y^^^^^-^ 



3 , -TPS xPT 



F 1 + 



1 



87r2F2 



BttiL 



(3.6) 



To conclude this section we integrate out the singlet field. In the infrared limit, ml ~ 
^ m|, XPT5 has to reduce to xPT, where the only dynamical degrees of freedom are 
the pions [2]. In oder to do so we keep m| fixed and expand the above observables around 
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(a) (b) (c) 

Figure 3: Diagrams contributing to the scalar field self-energy. Solid lines denote the scalar 
and the dotted ones pions. 



0. At NLO order in this expansion we indeed recover (3.6) after the identification 



h = - 2^73 - 3211 In 



■> 2 

rria 



1 



h =0^474 - 16r2 In -r^ C2m + C 



ml 



.2 

-Id 



(777,2 
mi 



(3.7) 



4. The scalar field two-point function 

In order to calculate the scalar field two-point function at NLO, we need to enlarge our 
initial set of operators in (2.19) to provide the needed counter-terms to renormalize it 

6C =f2pDSnS + d2md^Sd''S{x^U + xU^) + b2m.S'{x^U + xU^f + a2mS^{xh) 
+ e2^^'Re[det(x)]. 

Notice that the first of these counter-terms may be eliminated by using the LO equation of 
motion. We keep it for convenience, it renormalizes contributions to a contact term, to the 
scalar wave function, and to the scalar mass in the two-point function. The divergent parts 
of the low-energy constants in (4.1) are determined by the cancellation of divergences in the 
two-point function. The combination dim = 32&2m + 4a2m + 2e2m is indistinguishable in this 
quantity, and is renormalized as a whole. The divergent parts read as follows 

/2p:=/2p + r/A, dim:=dl^ + AA,{t = l,2), (4.2) 

with 

r/ = ^, Ai = 1^ (C2. - C2. + 6cL) , A2 = -||, (4.3) 

and A is defined as in (2.20). The set of diagrams contributing at NLO to the two-point 
function is shown in Fig. 3. The contributions at NLO to the off-mass-shell two-point 
function can be disentangled as 

^2pNLO = Us + Ps , (4.4) 
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(ms/ms.p,,ys,r/rpi,ys) 




Figure 4: Plot of the scalar mass (solid line) up to NLO and the sigma decay width (dashed 
line) at LO both normalized to their respective physical values as a function of the quark 
mass. and are set to zero, cf^ is taken from (4.10), and C2d and C2m from (5.8). 



where, as in pion two-point function, we find a contribution related to the unitarity cut. Us, 
and a polynomial term in s including the logarithm term Pg. Explicitly 



Us 



6cL 

F2 



J {ml, ml; s) {s - 2m^)' 



mis 



2 4 

S — 771 — 

where we have used the scale independent quantities dim 



(4.5) 



r 

2p 



327r2 



mz 



1 ^im 



A,, 
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(4.6) 



The scalar mass at NLO, defined as the pole of the scalar field two-point function, reads as 
follows 



ml 



m 



S, NLO 



327r: 



m 



327r 



2 2 

mtm% . - 



47r2 



2"2m 



6cL 
F2 



2\2 



(4.7) 



J(m2,m2;m|) (m| - 2m2) 



Contrariwise to [28] we do not find branching points in the behavior of (4.7) when the pion 
mass increases. In [28] two conjugate poles in s of the amplitude are found. The branching 
is a result of the different behaviour of these poles when they reach the real axis. In XPT5 
the denominator in the amplitude has a second order polynomial in s like in the lAM. The 
non-analytic piece has the same structure as the s-channel contribution to the lAM, but 
the t- and u-channel contributions are absend. They are treated as perturbations in xPT^. 
If we look for poles we find two complex conjugate poles like in [28]. However when we 
take into account the counting, we observe that the self-energy needs to be resummed when 
s — m| ~ m^/A^ only. As a consequence, only one pole remains, now with an imaginary 
part, related to the sigma decay width. 

If we set to zero the finite part of all the unknown counter-terms (/2p = d\^ = 0), the 
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scalar mass (4.7) scales almost linearly with the pion mass up to the threshold region, see 
full curve in Fig. 4, where the vt-tt production channel closes. 

The scalar decay width can be read from the unitary part U 5 once we plug the logarithm 
dependence of the J function and set it to be on-mass-shell 

r 3c?^ / Ami / 9 9 N 2 , . 

- = ¥—\'^ ^(ml-2ml) . 4.8 

2 87iF^ms\l ml^ ^ ^' ^ ^ 

Notice that it only depends on a single unknown LEG, ci^. Using the standard values for 
-F ~ and m.^, and taking specific values for the mass and width of the sigma resonance 
from [9] 

= 92.419 MeV , m^ = 139.57 MeV , mg, ccl = 441.2 MeV , Tccl/S = 272 MeV , 

(4.9) 

we obtain from (4.8) 

cL = 0.457. (4.10) 



5. Matching with lattice data: the pion mass and decay 

constant 

The expressions for the pion mass and decay constant (3.1) depend on several LEG not 
constrained by chiral symmetry, £3 , £4 , Cid , , C2m in addition to the quark masses and the 
bare parameters F, B, m|. At this point, and for fitting purposes, the finite part of the 
counterterms can be absorbed into F, B. Then, we have eight independent parameters 
at our disposal at NLO. 

Lattice QGD offers a new arena for determining the LEG. Unlike physical experiments, 
lattice calculations use different unphysical quark masses, providing for each point what can 
be considered as an uncorrelated experimental datum with Gaussian errors. We will use the 
lattice data based on maximally Uf = 2 twisted fermions to fit the LEG [13]. More precisely 
the data ensembles labeled as Al- A4, Bl- B6, Gl- G4, Dl and D2 in the Appendix G of 
that reference. Both finite volume effects and discretization errors are small in the data sets 
we use, and will be ignored in the following^. We will also use in our analysis a single lattice 
scale To = 0.446 fm for all data sets as a simplification, which is justified because its value 
varies very little from one data set to another. 

Given the limited quantity and quality of the available data, the number of free param- 
eters is too large to expect a brute force best fit to provide sensible values for all of them. 
We have rather used a general three-fold strategy: 

^This also implies that the isospin violation due to the twisted mass term in the lattice lagrangian is 
small, as it vanishes in the continuum limit, see [29]. Lattice data in ref.[13] corresponds to charged pions. 
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1. We identified Fps (3.1) with its physical value at the physical pion mass (4.9). This 
allows to write F as a function of the remaining parameters. This determination is 
done perturbatively. 

We have used the same procedure to fix the bare scalar mass through (2.14): one 
imposes to the tree level scalar mass to take its physical value (3.1). 

2. The parameter appears in our expressions always squared, hence we will use cf^ 
as the free parameter. As we have seen its value can be extracted through the sigma 
decay width, (4.10). 

3. For each point in the C2x {x = d, m) parameter-space we fitted the lattice data minimiz- 
ing an augmented chi-square distribution that includes both observables in (3.1) [30]. 
The augmented chi-square distribution is defined as the sum of the chi-square func- 
tions for each observable together with a set of priors for each of the free parameters 
to be fitted 



where g stands either for mps or Fps and gi for the corresponding lattice data at the 
quark mass m^. Furthermore, Xi refer to the fitted parameters being their total 
number. At NLO A^ = 3 and Xi = {-8/(2600 MeV) ,£3 ,^4}. The prior information on 
the LEG is obtained from naive dimensional analysis. For instance, if the parameter Xk 
is of order 0{1), we expect it to be in the range 0.1 < \\xk\\ < 10, which translates to 
setting ln(||afc||) = and ln(i?fc) = 1 for the kth parameter. We have taken logarithms 
in the prior functions to achieve equal weights for the subranges 0.1 < < 1 and 
1 < ||a:fc|| < 10. 

Within the previous outlined procedure we have fitted the expressions for the observ- 
ables mpg and Fps, (3.1), to the full range of available quark masses nig, the MS running 
mass at the scale /x = 2GeV. Throughout this section we have used A = 770 MeV as the 
renormalization scale value. The chi-square per degree of freedom, Xdo f is 



where A^' is the number of free parameters, including the fitted and scanned ones. Thus for 
the NLO fit A^' = 5. 




(5.1) 



Xd.o.f 




PS 



(5.2) 



2n - 1 - A^' 
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5.1. RESULTS 

Before we move ahead with the xP^^s case, and for comparison purposes, we present the 
outputs we obtain for xPT. From the previous steps in the fitting procedure we only make 
use of the first to fix the physical point of Fps and the third to introduce the priors to the 
corresponding parameters. 

5.1.1. Born approximation 

At LO, the bare parameter F is fixed by the physical decay constant F = 92.419 MeV and 
only S is a free parameter, obtaining 

5 = 2250.4 MeV, (5.3) 

„,;^-V, ^2 _ 1560 

With Xd.o.f - 

5.1.2. Next-to-leading order results 

For xPT the NLO expressions in (3.1) have three free parameters, B, Zg, and while F has 
been fixed perturbatively at the physical point. The fitting procedure leads to 

5 = 2499(10) MeV, = 0.91(6) x 10-^ = 7.13(5) x 10-^ (5.4) 

Using these values on the constraints imposed in the fist step of the fitting procedure we 
obtain F 

F = 86.36(1) MeV, (5.5) 

with the final value Xd.o.f — • 

The results, together with the lattice data, are plotted in dashed lines in Fig. 6. The 
adjustment for the pion mass to the lattice points is quite remarkable, with only a small 
deviation for large values of the quark masses. The pion decay constant fit also reaches a 
good agreement with the lattice data. 

The results for the LEG, (5.4), are compatible with standard values in the literature 

[33] [3 = 3.2 ±0.8, [34] [4 = 4.4 ±0.2. 

eq (5.4) [3 = 3.99 ± 0.04 , eq (5.4) [4 = 4.54 ± 0.01 . ^ ' ^ 

These estimates are in reasonable agreement with those obtained by resonance saturation 
[15] 

/3 = 2.9 ±2.4, /4 = 4.3 ±0.9. (5.7) 

The determination of the uncertainties in this section has been performed as follows. 
We assume each data point corresponds to Gaussian distribution with expected value and 
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variance defined by the data point value and uncertainty respectively, then we generate 
random data sets according to these distributions and perform a fit for each one. The final 
parameters are obtained from the average of the results of these fits, while the uncertainty 
is obtained from the variance. Comparing our results with those from Table 1 in [13], 



B = 2638(149)(132)MeV, F = 85.91(7)(t^^), k = 3.50(9)(;^o). h = 4.66(4)(;^3), we observe 



that F, and B values are within one sigma while I3 is within two sigmas. Note that our 
uncertainty analysis does not include systematic uncertainties because we have used only one 
set of data ensembles and we do not include finite size corrections. Statistical uncertainties 
in our fit are significantly smaller than those of [13]. This is because we have taken tq as a 
fix value rather than as an additional free parameter. If we estimate the uncertainty of tq 
as the one given in [13] and we extrapolate the effect to our results we obtain uncertainties 
in the same range as in [13]. Furthermore, since mpg depends quadratically on while Fps 
only linearly, our F and Z4 should be in better agreement with those of [13] than our B and 
I3, as it is the case. 

The estimation of uncertainties above is not directly applicable to the following section 
because for the fit to x PT5 expressions some of the parameters will be obtained by scanning 
a suitable range. In any case, we are not interested at this point in an accurate determination 
of the parameters but rather in finding out if parameter sets of this theory exists which 

are both compatible with lattice data and with physical observables. 

5.2. xPTs RESULTS 

The LO expressions are identical to those of standard xPT, therefore the same analysis 

as in the previous section applies. At NLO appear four extra free parameters to fit C2d , C2m , ^3 
and £4 and the non-analytical dependence on the light quark masses is greatly augmented 



The relative large number of free parameters appearing at NLO are an indication that 
there is no unique solution for the best fit. Indeed, if we look at the contour level plot of the 
Xd.o.f corresponding to the {c2d, C2m) region scanned, shown in Fig. 5, we can see regions of 
parameter sets with Xdof smaller than one. Thus any parameter set on those regions has 
to be considered a valid solution. Keeping this in mind, the following are the results for the 
best fit obtained, which have been used for Fig. 6 



(3.1). 



B = 1680.5 MeV , C2d = 1-21 , C2„ 
il = -1.12 X IQ-^ , £^ = 6.94x10 



-0.083, 



(5.8) 



,-3 



Using these in (2.14) and (3.1) we obtain the values of the remaining parameters 



F = 101.2 MeV, ms = 426MeV, 



(5.9) 
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Figure 5: xlo./ swept over a (c2d, C2m) grid corresponding to fits to NLO order expressions. 
The fits are forced to reproduce the pion decay constant, the mass of the sigma resonance 
and its width at the physical point. 
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Figure 6: The best fits of the LO (dashed line), NLO xPT (black solid line) and xPT^ (red 
solid line) expressions. Note that the LO expression is the same for xPT and XPT5. 
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rh (MeV) 


T (ci-':\i-\fc^rci-\ 
IjcLL L. I bLcLL )l bybl I 




Beyond NLO xPT 


3.54(19)(17) 


XPT, LO fit 


4.32 


XPT, NLO fit 


4.39 


xPTs, NLO fit 


3.26 



Table 1: The values obtained for the light quark masses from our fits (bottom panel), at LO 
xPT and xPT^ are identical. In the central panel we show the lattice QCD results from 
Ref. [13] that use chiral extrapolations beyond NLO xPT . In the upper panel we show the 
values from lattice QCD at the physical point [31]. 



with 



16.7 



(5.10) 



5.2.1. Scalar resonance contribution to xPT low-energy constants 

It is instructive to show how the LEG £3 and l\ of XPT5 above compare with the stan- 
dard LEG of xPT. This is done through the matching formula (3.7). We obtain that the 
corresponding values of [3 and read 

[3 = 3.40, [4 = 5.16, (5.11) 

thus, /s is compatible with the literature values in (5.6) and ^4 is somewhat higher. From 
(3.7) we can easily find out the fraction of Z3 and Z4 that is exclusively due to the light scalar 
field by setting £3 and £4 to zero. It amounts to a 43% for [3 (with opposite sign) and to a 20% 
for Z4. This suggest that the impact of the singlet field in both /3 and /4 is quite substantial. 
Note that the contributions of the scalar field to these LEG comes entirely through loops, 
and hence have nothing to do with the tree-level contributions obtained in [15]. 



5.2.2. Quark mass determination 

The last application we have explored is the determination of the light quark masses, and the 
comparison with the results obtained from xPT and lattice QGD. Given a set of parameters 
the expressions for mpg (3.1) and (3.6), become a function of rh. Setting m|g to the physical 
value of the pion mass we can solve the equation to obtain the value of m at the physical 
point. The results obtained for rh for the best xPTs and xPT fits are displayed in Table 
1. The expressions used for light quark masses match the order at which the fit has been 
performed, at NLO the equation has been solved perturbatively. 
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6. S-WAVE TT-TT SCATTERING LENGTHS 



V/ 



(a) (b) 

Figure 7: Diagrams contributing to the scattering lengths. Sohd hues denote the scalar and 
the dotted ones pions. 

Let us next consider tt-tt scattering. The diagrams contributing to the scattering ampli- 
tudes are depicted in Fig. 7. Due to the presence of a novel contribution coming with 
a scalar particle in the intermediate state, we expect a LO correction to the xPT re- 
sults. In fact this new contribution allows to test the quark mass dependence of gs-E-K-, the 
scalar-pion-pion coupling constant, as outlined in [32]. Following this reference, we define 
g\^^ := — IGvr limg_^^2 (s — m|)too ■, being too the isospin zero S-wave amplitude. If we take 
the ratio 



gs-n 



„phys 



mo — 2mi 



''S ^'''TrJphys 

and increase the pion mass we find a smooth decreasing function that vanishes around 
TTiyr ~ 368MeV. While this number roughly matches fit D of Figure 7 of [32] , we do not find 
the strong quark mass dependence in gs-Kn displayed in that figure. 

Let us now turn to the evaluation of the scattering lengths. Their explicit expressions at 
LO are given by 

° 7rF2 V 32 2 4ml - ^1 ^1 

As we have already shown there is a new contribution coming from the scalar exchange. 
Notice that although the size of the 4m^ — m| denominator is similar to m| at the physical 
point, for certain values of the quark masses this quantity may become small, and hence the 
self-energy corrections calculated in section 4 should be included. However, this will not be 
necessary for the values of the quark masses considered in this paper. In the sequel we will 
elucidate the precise role of the scalar in the scattering lengths. 
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6.1. First estimates 



Using (4.9) and (4.10) into (6.2) we can compute the values of the scattering lengths. The 
results, displayed in Table 2, show that, for a^, x^^^s overshoots the experimental value 
by roughly the same amount as LO xPT undershoots it, whereas for xPTg is roughly a 
factor of three off the experimental value, namely much worse than LO xPT, which provides 
a number pretty close to it already at LO. 

This missmatch may be understood as follows. In the decoupling limit (m| ^ m^, p"^) 
the contribution in Fig. 7(b) gives 

Fig.7 (b) = -lclF\D,UD^U^) _^ \ {D.UD'^U^) ^ ^I^JD.UD^U^)' , (6.3) 

i.e. it reduces to a contact term which is proportional to /i in xPT. By direct identification 
one finds the value of the xPT constant in terms of the xPTs parameters 

/^o = 192^^^ . (6.4) 

Note that the usual Air suppression factors coming from loop integrals are absent in the tree 
level calculation above. It is easy to check that the last operator in (6.3) reproduces the 
scattering lengths (6.2) in the decoupling limit 

~ nF^ [32 + 2mf'y ' ^° " vrF^ ^16 ' ^^'^^ 

Using (4.9) leads to ~ 38, roughly 20 times bigger and with opposite sign than the 
standard NLO value for this quantity in xPT, h ~ —1.8 [34]. This indicates that a large 
negative value is expected for ii, and, consequently, that NLO contributions are going to be 
large, at least the ones related to the ii operator. 



6.2. Matching with lattice data 

The available lattice results for the S-wave scattering lengths use relatively large pion masses, 
which makes chiral extrapolations less reliable. In fact, until recently only calculations of Oq 
were available [36, 37, 38, 39, 40, 41, 42], and the only existing calculation of both Oq and Qq 
neglects the disconnected contributions to the latter [14]. Nevertheless we shall use lattice 
data of the last reference in order to get a feeling on how xPTs performs with respect to 
the S-wave scattering lengths. 

As we discussed in section 6.1, the S-wave scattering lengths of xPT^ at LO are fixed 
once we input the mass and the width of the sigma resonance in addition to the pion mass 
and decay constant. Their evolution with the light quark masses is given by that of the pion 
mass and the LEG C2m- By making a combined fit to Oq and we obtain the dashed red 
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rL/xp . ( st at ) ( sy st ) 




n r\AOQf A A \ {Oii\ 

— u.U4zy (^44 ) [^zo ) 


Beyond NLO xPT 


0.220 ± 0.005 


—0.0444 ± 0.0010 


XPT, LO 


0.159 


-0.0454 


XPT, NLO 


0.228 


-0.0405 


xPTs, LO 


0.275 


-0.0121 


XPT5, LO+£i 


0.210 


-0.0296 


Linear sigma model 


0.696 


-0.0404 



Table 2: Values obtained for the scattering lengths from xPT, XPT5 and the linear sigma 
model (bottom panel). LO expressions are fixed by the pion mass and decay constant for 
xPT, plus the sigma resonance mass for the linear sigma model, plus the sigma resonance 
width for xPT^. NLO xPT values are obtained by fitting li + 2I2 and I3 to the lattice data 
of ref. [14]. LO+£i xPTs is obtained by fitting ii and C2m to the same lattice data. In the 
central panel we show theoretical results that go beyond NLO x PT from ref. [34]. In the 
upper panel we show the values of the scattering lengths extracted from experimental data 
[35]. 

line in Fig. 8. We observed that for xPT^ provides a better description of data than LO 
xPT (dashed black line), but for a much worse one. As argued in section 6.1, large NLO 
corrections due to £1 are expected. We may estimated them by just adding its contribution 
to LO expression. If we fit ii, we obtain the dashed red line in Fig. 8, and the following 
numbers 

a° = 0.210 , al = -0.0296 , C2m = -0.443 , ii = 967i^ii = -16.9 . (6.6) 

Note that we get a large negative number for ii, consistent with the expectations. Notice 
also that the value of C2m above justifies the use of formula (6.2) for aj] (i.e. with no self- 
energy corrections in the m| — 4m^ denominator), with the possible exception of the point 
corresponding to the largest pion mass. We see that the description of both scattering lengths 
improves considerably, the quality of a^ being comparable to that of NLO xPT (black solid 
line). The plots of NLO xPT in Fig. 8 are obtained by fitting l[ + I2 and /g. The values 
delivered by the fit are 

ri + 2r2 = 1.4 , Is = -9.3, (6.7) 

which differ quite a lot from the standard values in xPT at one loop, for instance, I1+2I2 ~ 9.0 
is given in [34] and I3 ~ 2.9 in [2]. In fact if I3 is fixed to the last value rather than fitted a 
very bad description of aj] is obtained, whereas the one of al remains quite good. 

The results above encourage us to attempt an extraction of the sigma resonance param- 
eters from the lattice data. We obtain from the fit (to both al and Og) 

C2,„ = -0.228 , £1 = -10.9 , c?rf = 0.304 , = 483MeV , (6.8) 
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Figure 8: The best fits of the LO xPT (black dashed hne), NLO xPT (black solid line), LO 
xPT^s (red dashed line) and LO x^T^s augmented by the operator proportional to ii (red 
solid line). Red dots are lattice data from [14]. 

which produce the following numbers for the sigma decay width and the S-wave scattering 
lengths 

nis = 486MeV , - = 236MeV , a° = 0.177 , = -0.0361 . (6.9) 

The numbers above are quite reasonable for a LO approximation augmented by ii, even more 
if one takes into account that the lattice data is at relatively large pion masses. It shows that 
our approach may eventually allow for a precise extraction of the sigma resonance parameters 
from lattice QCD. Note in particular that the value of C2m is compatible with the region of 
low Xdo f of PiS- 5 and that ii remains with a large negative value. 



7. Discussion and Conclusions 

We have considered the possibility that the spectrum of QCD in the chiral limit contains 
an isosinglet scalar with a mass much lower than the typical hadronic scale A^, and have 
constructed the corresponding effective theory that includes it together with the standard 
pseudo-Goldstone bosons, xP'~^s- This effective theory has the same degrees of freedom 
as the linear sigma model, but differs from it in two important points. First of all, it is 
conceptually different because the mechanism of spontaneous symmetry breaking is assumed 
to occur at the scale A^, and hence it is not described within the effective theory. Second, 
there is a power counting and hence the LO Lagrangian can be augmented at the desired order 
by adding power suppressed operators. The LO Lagrangian has initially four free parameters 
more than the linear sigma model, and hence enjoys a larger flexibility to describe data. As 
explained in the section 2.3, one of these parameters (cim) must be set to zero for consistency, 
whereas in the linear sigma model it takes a non-zero value. If we force the LO fits to the 
pion mass and decay constant to go through the linear sigma model values we obtain a 
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Xdof^ 135, namely worse than in LO xPT^ (which coincides with LO xPT). Inputing the 
sigma mass in [9], the hnear sigma model delivers a relatively low value for the decay width 
(r/2 = 188), a very large value for the isospin zero scattering length (cq = 0.696) but a 
pretty reasonable one for the isospin two one (cq = —0.0404), see Table 2. 

At tree level xPT^ gives definite predictions for S-wave scattering lengths if the mass 
and decay width of the sigma resonance are used as an input, which are shown in Table 2. 
Neither the value of the isospin zero one (og) nor the one of the isospin two (oq) are close 
to the experimental numbers. Although the value of is slightly closer to it than the one 
obtained in tree-level xPT, the value of is much further away. As argued in section 6.1, 
this is due to the fact that sizable NLO corrections due to a large value of ii are expected. 
If we simulate them by letting ii be a free parameter, the combined fits to the lattice data 
of ref. [14] to Qq and Oq become rather good, see Fig. 8. Note that, although NLO xPT 
produces a better description of and Oq if + 211 are fitted to data, the values 

delivered by the fits of those LEGS are incompatible with the ones currently used in xPT. 
We have also shown how the combined fits to the S-wave scattering lengths may be used 
to extract the resonance parameters of the sigma from chiral extrapolations of lattice QCD 
data. 

Loop corrections in xPTs have been explored in the calculation of -F^r and m,r at NLO. 
The dynamical scalar field introduces new non-analyticities in the quark mass dependence 
of these observables, and requires a renormalization of B and F, which are absent in xPT. 
The fits to the lattice data of ref. [13] for these observables at NLO in xPTs are of similar 
quality as those at NLO in xPT. However, when the value of the average light quark masses 
is extracted from the fit, xPT^ produces numbers that are closer to those of direct lattice 
extractions than xPT does, see Table 1. The self-energy of the scalar field has also been 
calculated at NLO. 

We have restricted ourselves to the fiavor SU{2) case, the extension to flavor SU{3) 
is straightforward. In fact because flavor is conserved at any vertex, the contribution to 
observables with pions involving scalar fields in internal lines are identical and independent 
of the group, at the order we have calculated. Furthermore, because we will have more 
parameters at our disposal and ms ~ niK ~ we expect that the tension between the 
different contributions to higher chiral orders [43] is alleviated. 

Let us also mention that Lagrangians identical to the first line of (2.8) are currently being 
used in the context of composite Higgs models [44]. In that context, XPT5 would correspond 
to an effective theory at the electroweak scale under the assumption that the spontaneous 
symmetry breaking mechanism takes place at a much higher scale. Small explicit breaking 
of custodial symmetry at that scale may be taken into account by terms similar to those in 
the second line of (2.8). 

In summary, we have shown how to consistently introduce a light isosinglet scalar particle 
in a chiral effective field theory framework, XPT5. This has consequences concerning the 
dependence of physical observables on the light quark masses, which have been shown to 
be compatible with current lattice data. We have also shown that our formalism has the 
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potential to extract the mass and width of the sigma resonance from lattice QCD data. 
Finally, it would be interesting to explore the consequences of x^T^s in the chiral approach 
to nuclear forces [45] (see [46] for a recent review), since the exchange of a scalar particle 
is known to be an important ingredient of the nuclear force in one-boson exchange models 
[47]. 
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8. Appendix 



Through the calculations we have used the following set of integrals 

A^-*^ r df^k 1 



i J (27r)°' fc2 — + ie 



K (-2A + lla) . 



1 



I J (27r)°' — + ie [k — pY — ml + ie 
+ J[ml,ml;p^] . 



Which finite parts are given in terms of 
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(/ + uY - A=^ 
(p2 - uY - A2 



(8.4) 



with A = Aab = ml — ml, S = ml + ml, and u'^ = {p^ — {ma — mh)'^) {jp' — {ma + m^Y)- 
The expression in 8.4 is correct in the momentum region p^ < {ma — mYY- The analytic 
continuation to higher momentum regions is obtained using the following prescription p^ — > 
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+ ie. As a last comment, we have used MS subtraction scheme 



(8.5) 
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